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Abstract 

The study of Coulomb excitation of 11 Be to the first excited state in intermediate energy colli- 
sions with heavy targets is presented. The existing experimental data are reanalised by including 
the probability of projectile survival into the calculation of Coulomb excitation cross sections. 
The survival probabilities are calculated using a recently developed global optical model potential 
tailored in line with the double folding model. The extracted B(E\) values for the transition 
^ + — ► ^ in n Be are found to be slightly larger than those obtained so far using the fro-recipe in 
Coulomb excitation theory. 
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The Coulomb excitation process has been proven to be a useful tool to study nuclear 
structure (see, e.g., Ref. l]]). It is an attractive test of the different nuclear models as the 
Coulomb interaction is very well known. Recently a number of experiments has been done 
to study the structure of light neutron rich nuclei, such as n Be |2|, yl, |4( , 15 C, 17 C, 19 C and 
several oxygen isotopes, using the Coulomb excitation process j^. 

The isotope n Be is one of the most studied halo nuclei nowadays. The reason why 11 Be 
is an attractive system to study by Coulomb excitation lies in the peculiar structure of the 
spectrum of this nucleus. Both the | + ground state and the only bound excited state | 
are weakly bound and have energy difference of 0.32 MeV, thus making the El transition 
in 11 Be the fastest known between bound states. 

There have been several experiments that measured the excitation of the first excited state 
in n Be and extracted the corresponding B(E1) value. The adopted value ca 0.116 e 2 fm 2 
was obtained by Millener et al. in Ref.p] by averaging the results of three experiments on 
the lifetime of the excited state using a Doppler-shift technique. 

The experiment done in GANIL [2] studied the inelastic scattering of 11 Be on a lead 
target at a bombarding energy of 45 MeV/u. The obtained cross section was only 40% 
of the one predicted by calculations for pure Coulomb excitation. To explain this large 
difference it was proposed that the higher order effects may contribute, but calculations 
done in Ref.Q showed that the cross section falls by only 4% if coupling to continuum was 
taken into account. The inclusion of monopole and quadrupole modes of nuclear excitation 
also performed in the work ?J resulted in the slight increase of the cross section less than 
2%. 

Similar experiment was done by Nakamura et al. in Ref. [3] at E = 64 MeV/u. The B(E1) 
value extracted from the cross section assuming pure Coulomb excitation was comparable 
with the result of the Millener analysis. Later the inelastic scattering of 11 Be was measured 
at MSU y on lead and gold targets for energies 60 and 58 MeV/u, respectively. The 
extracted B(E1) value confirmed that of [^| and agreed, at least marginally, with the lifetime 
experiments. For summary of this see Table |Hj 

In the analysis performed in the mentioned papers, the authors used the formalism of pure 
Coulomb excitation and excluded nuclear processes approximately by using a low impact 
parameter cutoff bo . In this work we show that the inclusion of the full fledged survival 
probability ^(fr)) 2 in conjunction with b , though slightly improves the agreement with the 
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data, it does allow one to perform a more realistic smooth cutoff based analysis. 

In the calculation of the Coulomb excitation cross section, the important ingredient of 
the model is the minimum value of the impact parameter, 6 > from which the integration 
of the excitation amplitude is performed. This b Q is roughly determined by the sum of the 
projectile and target radii. n Be, being a halo nucleus, has a very diffuse structure, which 
makes the definition of the minimum impact parameter obscure. 

To improve upon the calculation we multiply the probability of the Coulomb excitation at 
a given impact parameter by the survival probability ^(fe)! 2 calculated for the system under 
consideration. In this work we take the S-matrix from the optical model. The calculation of 
the elastic S-matrices requires the knowledge of the densities of the projectile and the target. 
The halo density of n Be was calculated in the framework of the particle-rotor model, which 
includes the excitation of the rotational 2 + state of the core 10 Be 

As was mentioned above, the contribution of nuclear excitation is small. This is so 
since nuclear effects are limited to a very small impact parameter region around the grazing 
value. In fact, as shown in Refs. the Coulomb excitation was found to be by far 



the dominant piece of the cross section. Further, in Ref.[4j an experiment was done with 
light targets to study the importance of nuclear excitation. The results were 4.0 mb for 
carbon target and 1.7 mb for beryllium target. Since the nuclear cross section corresponds 
to the area of a ring around grazing b , it scales with A as A 1 / 3 . Accordingly, we find for the 
lead target approximately 10 mb nuclear excitation contribution implying a mere 3% effect. 
Therefore in the following we ignore the nuclear excitation effect. 

The semiclassical model is usually used to describe Coulomb excitation at intermediate 
and high energies. This model assumes the straight line trajectory for the projectile and 
treats quantum mechanically the absorption of radiation by the nucleus. The formalism 
of this method was presented in Ref.ll2« and later was extended to the relativistic case in 

n 

Ref.HJ. For the hieh energies the first order perturbation theory is a good approximation 
to calculate the amplitudes for Coulomb excitation. In the first order perturbation theory 
the process of Coulomb excitation can be described as emission and absorption of virtual 
photons [14|. We included only El multipolarity in our analysis since for the transition 
| + — > | in n Be the dipole multipolarity is the dominant one Q|. 

Using the formalism of virtual photons the Coulomb excitation cross section has the 
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following form: 



a c = f° ^dVl = 2tt r P{b)bdb (1) 
Jo di I J bo 

with 

P(b) = ^B(El)N E1 (u : b), (2) 

and bo = acoty, where a is half the distance of closest approach for head-on Coulomb 
collision. Nei{uj, b) is the number of virtual photons given by 



_ j , Zfa ( A /c\2 



Kl (0 + ^ (0 



(3) 



7T 2 V 6 / Vt>, 

where Z x e is the charge of the target, £ = to = E ex /h, E ex is the excitation energy of 
the state in the projectile, v is the relative energy, 7 is the relativistic factor and a is the 
fine structure constant. The behavior of the integrand in the Eq.Q is determined by the 
impact parameter dependence of the modified Bessel functions of the zero and first order 
K and K\. 

Coulomb recoil was taken into account using the method of Winther and Alder 
replacing b in the expression for P(b) by b' = b+7ra/2j. Further, we multiply P{b) of Eq.JU) 
by the survival probability ^(fo)! 2 , calculated from the optical potential NLM3Y of Ref.|16| 
and then write 

P(b) = P(b).\S(b)\ 2 . (4) 

The optical potential was calculated from the double folding model, using projectile and 
target densities and the effective nucleon-nucleon interaction. The effect of the exchange- 
related non- locality is taken fully into account. The local-equivalent, energy-dependent 
potential can be written as |16] 

V(r, E) = V N (r, E) - iW{r, E), 

where 

V N (r,E) = V F (r)e-^ 2 / c \ 
where c is the speed of light and v is the local relative velocity between the two nuclei 

v 2 (r, E) = -(E- V c (r) - V N (r, E)). (5) 
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The local folding potential Vp(r) is obtained in the usual way. For the Coulomb part, 
Vc(r) we employ a similar folding prescription using the proton densities. The imaginary 
part of the optical potential was found to be 

W(r,E)&0.78V N (r,E). (6) 



The determination of the number 0.78 is discussed in Ref.[17|. This number was obtained by 
fitting calculated elastic scattering cross sections with experimental data for the following 
systems: 12 C+ 12 C, 16 O, 40 Ca, 90 Zr, 208 Pb; 16 O+ 208 Pb, 40 Ar+ 208 Pb and for the wide range of en- 
ergies. The calculations showed that the number 0.78 is approximately system-independent. 

We employ the same potential here by using the appropriate halo+core density of n Be. 
For the densities of the heavy targets and the two-parameter Fermi distributions were used 

P(r) = - P ° r _ Ro , (7) 
1 + exp ■ L -^ M 

where a = 0.56 fm, R = 1.3L4 1 / 3 - 0.84 fm and on is determined by the normalization 



condition. These densities were obtained in Ref.|16j] with the aim of providing a global 
description of the nuclear interaction, based on an extensive study involving charged distri- 
butions extracted from electron scattering experiments and theoretical densities calculated 
through the Dirac-Hartree-Bogoliubov model. For the density of the core of n Be, 10 Be, we 
used the Eq.® with a = 0.5 fm and R = 1.75 fm, which gives i? rms ( 10 Be) = 2.3 fm. 

We use the existing particle + core excitation model to describe the halo structure of a 
light halo nucleus such as 11 Be jsj]. The inclusion of the core excitation allows for the coupling 
between the collective degrees of freedom of the core and the orbital motion of the single 
neutron. In this calculation we assumed a quadrupole- deformed core (with deformation 
parameter fa) and included only the ground state + and the first excited rotational state 
2 + . The interaction between the valence particle and the core was described by a deformed 
Woods-Saxon potential. The spin-orbit interaction, V so , is the standard undeformed one. 
To reproduce the adopted B(E1) of n Be we adjusted the model parameters obtained in 
Ref. 8j and these parameters are presented in Table HI The particle-core model with these 
parameters gives a radius for n Be of 3.0 fm, the B(E1) = 0.116 e 2 fm 2 . 

In FigJT] we show the calculated squared moduli of the elastic S-matrices for 11 Be + 
208 Pb at 59.7 MeV/u. It is seen that using the density with the halo reduces the survival 
probability. 
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TABLE I: Parameters of the particle-core model used to calculate density of 11 Be. V so and 
Vw S VEN ^° DD are the strengths of spin-orbit and Woods-Saxon potentials, a and R a are the diffuse- 
ness and the radius of the potentials. is the quadrupole deformation parameter and E*{2 + ) is 
the excitation energy of the core. 

In Ref. jl^ the projectile survival probability was determined from the angle dependence 
of the Coulomb dissociation cross section. The very small number of experimental points (3 
in the grazing region) with their rather large error bars is certainly insufficient to make mean- 
ingful comparison. In Ref. [181] the survival probability was approximated by the function 
with the shape 1/(1 + exp (— (6 — bo) /a)) with parameters bo = 12.3(1.2) fm and a = 0.9(0.6) 
fm. In fact, one should include the size of the error bars in the presentation of the experimen- 
tal survival probability. For the purpose of completeness we show in FigElthe experimentally 
determined survival probability of Nakamura et al. presented in the shaded area together 
with ours shown as the full line. It is clear that our calculation of the survival probability 
agrees reasonably well with Nakamura's one. 

The experimental data on Coulomb excitation of 11 Be is summarized in Table ILT1 The fifth 
column shows the values of B(E1) extracted from the corresponding Coulomb excitation 
cross section by assuming pure Coulomb excitation and choosing appropriate fen, that is 
without including projectile survival probability. In the analysis of Nakamura et al. [3j Eq.(£Q) 
was used with b = 12.3 fm which was determined from the impact parameter dependence 
observed in the Coulomb breakup of n Be ^j]. Fauerbach et al. in Ref.J^J used the similar 
formalism to extract B(E1) values from the experimental cross sections with bo obtained 
from the aperture of the experimental setup. 

In our calculation we use P(b) in Eq.(J^) and start the integration using the experimental 
value of bo, which is determined by the aperture of experimental setup. Thus for the case 
of, e.g., n Be + 208 Pb at 59.7 MeV/u we integrated the probability shown by dotted line 
starting with b = 11.5 fm. 

The theoretical cross sections calculated in this work are presented in columns 6, 7 and 
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FIG. 1: The projectile survival probability \S(b) | 2 calculated for n Be + 208 Pb at 59.7 MeV/u. 
The solid line shows the projectile survival probability calculated using the densities of the halo, 
dashed line stands for the calculation using the global densities of the Eq.(0). 
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FIG. 2: The projectile survival probability \S(b)\ 2 calculated for n Be + 208 Pb at 72 MeV/u 
including halo density. The shaded area represents the projectile survival probability determined 



from experiment in 
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8 of Table |HJ First we assumed the B(E1) values from the previous analysis and 
calculated Coulomb excitation cross sections without taking into account Coulomb recoil and 
using density of Eq.(jJJ). We obtained the cross sections slightly larger than the experimental 
for n Be + 208 Pb at 64 and 45 MeV/u as we integrated the probability from 6 obtained from 
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FIG. 3: The impact parameter dependence of P(b)b and P(b)b. The case of n Be + sosp^ 
59.7 MeV/u is considered. P(b)b was calculated with the halo density (solid line) and with global 
densities of the work [16( (dashed line) in \S(b)\ 2 , see text for details. 

ri 

experimental aperture and not from grazing impact parameter as was done in |3[. Taking 
into account Coulomb recoil reduced the calculated cross section by approximately 2%. The 
further reduction (approximately 2%) in the calculated cross section was obtained by using 
the halo densities of n Be. 

In FigOHwe show the impact parameter dependence of the Coulomb excitation probabil- 
ities P(b) and P(b) for the case of n Be + 208 Pb at 59.4 MeV/u. The inclusion of projectile 
survival probability reduces the Coulomb excitation probability for impact parameters from 
11.5 fm and up to 17 fm. This results in the reduction of the cross section by 4% compared 
to the result of the pure Coulomb excitation with b = 11.5 fm. For the case of n Be + 197 Au 
the the reduction of the calculated cross section was 5%. 

Finally, we extracted the revised values of B(E1) using P(b) and taking Coulomb recoil 
into account. We found no difference between the values of B(E1) extracted in the work 
and results of our analysis. The revised B(E1) values obtained from the experimental cross 
sections of the work fjj are approximately 4% larger than those obtained in [J] . The revised 
B(E1) values are presented in the last column of the Table ITT1 

In conclusion, a better treatment of nuclear absorption exemplified by the use of an 
appropriate survival probability in the calculation of the Coulomb excitation cross section 
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Ref. Target Beam energy a c (exp.) B(E1) (old) a c (th. (1)) ct c (th. (2)) <r c (th. (3)) B(E1) (this work) 
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TABLE II: The experimental and theoretical Coulomb excitation cross sections for different en- 
ergies and extracted B{E1). The fifth column shows the B(E1) extracted from experimental a 
using pure Coulomb excitation in Refs.jj], 0|. Columns 6,7 and 8 show calculated cross section 
(using B(E1) from column 5) assuming no Coulomb recoil and no halo-type density (th. (1)); with 
Coulomb recoil, but no halo- type density (th. (2)); with Coulomb recoil and halo- type density 
(th. (3)). The last column shows the revised values of B(E1) which reproduce experimental cross 
section using P(b) and taking Coulomb recoil into account. 

of 11 Be leads to slightly smaller cross-section and then a slightly (about 4%) larger B(E1) 
value. We expect similar effect in the Coulomb dissociation cross-section of 11 Be and other 
halo nuclei. 
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